Efficient Quantum State Estimation by Continuous Weak Measurement and 

Dynamical Control 
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We demonstrate a fast, robust and non-destructive protocol for quantum state estimation based on 
continuous weak measurement in the presence of a controlled dynamical evolution. Our experiment 
uses optically probed atomic spins as a testbed, and successfully reconstructs a range of trial states 
with fidelities of ~ 90%. The procedure holds promise as a practical diagnostic tool for the study 
of complex quantum dynamics, the testing of quantum hardware, and as a starting point for new 
types of quantum feedback control. 
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Fast, accurate and robust quantum state estimation 
(tomography) is important for the study of complex 
quantum systems and dynamics, and promises to be an 
essential tool in the design and testing of hardware for 
quantum information processing . Previous demon- 
strations range from optical to atomic and molec- 
ular systems 4], but with few exceptions these proce- 
dures have proven too cumbersome to be of use as prac- 
tical laboratory tools. The procedure of quantum state 
estimation is usually formulated in terms of strong mea- 
surements of an informationally complete set of observ- 
ables. Each such measurement erases the original quan- 
tum state, so the ensemble must be reprepared and the 
measurement apparatus reconfigured at each step. Here 
we demonstrate a general approach based instead on con- 
tinuous weak measurement Using a weak measure- 
ment spreads quantum backaction across the ensemble 
and dilutes it to the point where it does not significantly 
affect any individual member. In the absence of back- 
action the quantum state remains largely intact, subject 
only to minimal damage from errors in the external drive 
fields and coupling to the environment. This allows us 
to estimate the state in a single interrogation of the en- 
semble, based on the measurement of a fixed observable 
O and a carefully designed system evolution. 

In the Heisenberg picture the chosen dynamics leads 
to a time-dependent observable, O — > 0(i), and the mea- 
surement history can be made informationally complete 
if the system is controllable, i. e. the dynamics can gen- 
erate any unitary in SU (d) where d is the dimensional- 
ity of Hilbert space. With a non-destructive measure- 
ment and near-reversible dynamics, the entire ensemble 
remains available at the end of the estimation procedure, 
in a known quantum state that can be restored close to 
its initial form if desired. In principle, this means that 
the knowledge gained can be used as a basis for further 
action, for example real-time feedback control Q or error 
correction [l| . From a practical viewpoint our procedure 
is highly efficient: the ~ 1ms interrogation time is lim- 
ited only by the control and measurement bandwidths, 



and data analysis is performed offline. It is also robust, 
in the sense that imperfections in the experiment can be 
included in the analysis if known, or estimated along with 
the state if they fluctuate in real time. 

The quantum system used in our laboratory imple- 
mentation is the total spin-angular momentum (electron 
plus nuclear) of an atom in the electronic ground state, 
here the F = 3 hyperfine manifold of 133 Cs. A detailed 
theoretical description of our estimation procedure and 
how to model this system can be found in 5]. Our Cs 
ensemble is prepared by laser cooling and optical pump- 
ing, and consists of a few million atoms in a cloud with 
radius ~ 0.3 mm. The spins are coupled to the polar- 
ization of an off-resonance probe beam that serves as a 
meter for the measurement and can be read out in an 
arbitrary basis with a shot-noise limited polarimeter Q , 
as illustrated in fig.^i. The spin-probe coupling depends 
on the atomic tensor polarizability, and can be expressed 
in terms of the single-atom spin F and the probe Stokes 
vector S, 
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where the coupling constants & depend on the probe 
frequency, the cloud optical density, and the atomic tran- 
sition Q. Terms containing Si, S2, and S3 lead to rota- 
tion of the Stokes vector around the 1, 2 and 3-axis of the 
Poincare sphere, by angles proportional to the ensemble- 
averaged values of the respective spin observables. We 
use a probe initially polarized along the x-axis, corre- 
sponding to an input Stokes vector S^ n = So (1,0,0), 
and measure the induced ellipticity, or S3 component of 
the output. For the moderate optical density clouds used 
here all angles of rotation are small, the polarization mea- 
surement is sensitive only to rotations around the 2-axis, 
and the measured observable is O — F x F y + F y F x . Cou- 
pling the entire spin ensemble to a common probe can, in 
principle, produce correlations between individual spins 
but this effect is negligible for our optical densities. 
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FIG. 1: Continuous measurement and quantum state estimation for atomic spins, (a) Schematic of our experiment. PBS: 
polarization beamsplitter; PD: photodetector; A / 4: quarter- wave plate. The control magnetic field has constant magnitude 
and is confined to the x-y plane; the time dependent angle <j> is optimized for best estimation fidelity, (b) Simulated (dark blue) 
and observed (light red) measurement signals for test states \rriF — —3) (top), \ip c ) (middle) and pmix (bottom), (c) Input and 
estimated density matrices (absolute values) corresponding to the simulated and observed measurement signals. (Color online) 



To extract information we evolve the spins to ob- 
tain a time varying observable O (t) . Coarse-graining 
over time yields a discrete set {Oi} and a measure- 
ment record {MJ, where =Tr[0; p ] + AM. Here 
po is the quantum state to be estimated, and AM is a 
Gaussian white noise process. Full system controllabil- 
ity is required to make {Oi} informationally complete, 
and can be achieved with a Hamiltonian of the form 
H s (t) = gf fiB B (t) • F + (3 Hj s F% produced by a time- 
dependent magnetic field and the probe-induced tensor 
light shift (eq. Note that for alkali atoms with 

an nSi/2 electronic ground state, the tensor light shift is 
proportional to the photon scattering rate 7 S and there- 
fore comes at the cost of decoherence. By probing Cs 
atoms on the Di line we achieve (3 = 8.2, enough to 
allow significant coherent evolution before optical pump- 
ing affects the spin state. We model our system dynam- 
ics with a master equation that fully accounts for opti- 
cal pumping and experimental imperfections. The latter 
include bandwidth limits for our magnetic coil drivers 
and detection system, a 6% spatial inhomogeneity in the 
magnitude of the light shift which can be independently 
estimated, and a ~ 1% drift in our magnetic coil calibra- 
tions which must be included as free parameters in the 
estimation procedure. 

We have evaluated the performance of our state es- 
timation procedure by applying it to a range of test 
states. Nearly pure test states with = 0, ±F were 
produced by optical pumping, and verified by Stern- 
Gerlach measurements of the magnetic populations. Ge- 
ometric rotations of these states were obtained by Lar- 
mor precession, and more general test states were pro- 
duced by evolving the system with a magnetic field and 



the tensor light shift. A strongly mixed state was ob- 
tained by omitting the optical pumping step and work- 
ing directly with the laser cooled sample. Figure QJd 
compares the measurement record observed in a sin- 
gle ensemble interrogation to the one predicted by our 
model, for three representative inputs corresponding to a 
spin-polarized state \ttif = —3), a coherent superposition 
\i/j c ) = A= (\rriF = 3) + i \rriF = — 3)), and a nearly max- 
imally mixed state pmix ~ I / 7. Two key features are 
immediately apparent: the measurement records differ 
substantially for different input states, and the observed 
and calculated measurement records agree very closely. 
Both conditions are clearly necessary for the procedure 
to generate unique and accurate state estimates. 

Based on our model we can determine the prob- 
ability distribution for the measurement record, 
P ({Mi} | {po, Pk}), conditioned on the initial quantum 
state po and any unknown experimental parameters 
Pk- We use this to obtain a maximum likelihood 
estimate pML for the state, corresponding to the peak 
of the posterior distribution, P({po, Pk} \ {Mi}) = 
AP({M Z } | {p , Pk}) P({p^Pk}), where P({p , Pk}) 
is the prior information and A is a normalization 
constant. Here we include as prior information only that 
po must be a physical state, i. e. a positive Hermitian 
operator with unit trace. The conditional distribution 
P ({Mi} | {po, Pk}) is a multi-parameter Gaussian 
centered on pLs, the ordinary least-squares fit to the 
measurement record. Its covariance matrix determines 
the uncertainty in each of the (2F+1) 2 real- valued 
parameters required to specify po, and consequently how 
well an arbitrary state can be estimated. In setting up 
the procedure, the control magnetic field was optimized 



to minimize these widths and reduce their sensitivity to 
errors. In practice the state estimate is obtained from 
the measurement record in a two-step process. First we 
determine the least-squares fit pls, which is frequently 
non-positive due to mismatch between the modeled 
and actual dynamics, and then we find the positive 
density matrix that lies closest to pLS- This is our 
maximum likelihood estimate Pml- Figured; compares 
density matrices for the input and maximum likelihood 
estimates for each of our three representative test states. 
Visual inspection clearly shows the excellent agreement 
achieved by our procedure. A more quantitative per- 
formance measure can be obtained by calculating the 

fidelity T = ^Try^ p~J 2 pML Pc/ [HI which for these 

three states fall in the range 0.89 - 0.98. Averaging over 
our entire sample of test states we achieve a mean fidelity 
T ~ 0.86 for estimates based on a single measurement 
record, and a slightly improved T « 0.91 for estimates 
based on the average of 128 measurement records. 

An alternative visualization of spin states can be given 
in terms of Wigner function representations . Figure 

shows the estimated states for six discrete time steps 
during Larmor precession of an initial state \mp = 3). 
As expected for these so-called spin-coherent states, the 
Wigner functions resemble Gaussian "wavepackets" mov- 
ing on a sphere. A slight but visible sheering of the quasi- 
probability distribution occurs due to the F% term in the 
probe-induced light shift, which is expected since this is a 
well known mechanism for spin squeezing 0. Figure 
shows estimated states at later times when the effects of 
squeezing are more pronounced. As the squeezing ellipse 
wraps around the sphere it eventually evolves into the 
coherent superposition \ij) c ) Q, the state used as input 
in fig. n^&c. This series of snapshots demonstrates that 
our technique can be used to observe and visualize quite 
complex dynamical evolution. 

With known input states we can observe how the state 
estimate converges as the measurement record builds up. 
Figure shows examples of the evolving fidelity and state 
estimate for pure and mixed inputs. In all cases, the es- 
timation starts with an unbiased guess, corresponding to 
a maximally mixed state. At short times, the dynam- 
ics has explored only a small part of the set {Oi}, and 
the measurement record provides good information only 
about a few of the parameters that go into specifying the 
state. As a result the estimate fluctuates randomly with 
each time step, subject to noise and model/experiment 
discrepancies. With increasing time, the dynamics ex- 
plores a larger set of observables, and the probability 
distributions are narrowed for every parameter. At ap- 
proximately 250 ps the constraints become strong enough 
that only a limited family of similar density matrices are 
consistent with the measurement record, at which point 
the estimate pml stabilizes and the fidelity improves. 
Figures |3Ji&b show the behaviors for inputs \mp = —3) 
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FIG. 2: Wigner function representations of an evolving spin 
state, (a) Snapshots of the quantum state estimate as an 
initial \ijif = —3) state undergoes Larmor precession in the 
presence of the probe. The nominal precession angle is in- 
dicated in each case, (b) Pronounced squeezing of the spin 
wavepacket occurs at later times, and at = 2868° the sys- 
tem has evolved into a coherent superposition \ip c ) of spin- 
up/down wavepackets. The viewpoint is fixed, except for the 
snapshots at = 630° where it has been changed to show 
the wavepacket on the opposite side of the sphere, and at 
6 — 2868° where it has been adjusted to show both lobes of 
\ip c )- (Color online) 



and |^ c ), where the transition at 250 ps brings the fi- 
delity close to its final value. The details of the jump 
depend on the semidefinite program solver used to en- 
force positivity, but it is always dramatic for pure states 
where the positivity requirement excludes most nearby 
points in parameter space. Fig. [3^ shows the behavior 
for an input pmix, where most nearby points in param- 
eter space fulfill the positivity requirement. As a result 
the initial jump is less dramatic and the fidelity increases 
more gradually to its maximum value at the end of the 
measurement period. 

In the Heisenberg picture the state po is time indepen- 
dent, and the estimate pML evolves only because more in- 
formation becomes available in the measurement record. 
Once the measurement is complete the estimate pML also 
becomes time independent. Because the dynamics are 
known we can then transform to the Schrodinger pic- 
ture, where the estimate becomes valid for any time, 
po (t) « pml (t). During the measurement period there 
is, of course, a gradual and irreversible — but entirely 
predictable — loss of information from the randomizing 
effects of optical pumping and drive field errors. As a 
result an initially pure state will become mixed, and the 
damage can be quantified by keeping track of the largest 
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FIG. 3: Evolving fidelity, state estimate and purity. Plots 
show the fidelity (light red) of the estimate pML and the 
largest eigenvalue (dark blue) of the evolving state po as 
a function of elapsed measurement time, for inputs (a) 
\rriF = —3), (b) \ip c ) and (c) p m ix- Inserts show the estimated 
density matrices (absolute values only) at a few representative 
times. (Color online) 

eigenvalue of po (t) . Figure 01 shows a small decline for 
the input ~ \rriF = —3), a slightly larger decline for the 
more fragile superposition \i/j c ), and essentially no dam- 
age to the mixed state which is already randomized. Ir- 
respective of the input state, optical pumping transfers 
~ 10% of the ensemble into the F = 4 manifold during 
the 1.5 ms interrogation, but these atoms are invisible to 
the probe and do not contribute to the estimate or its 
purity. 

With this work we have shown that fast, accurate 
and non-destructive quantum state estimation can be 
achieved with continuous weak measurement and dynam- 
ical control. In our current implementation the estima- 
tion fidelity is limited by imperfect control of the slowly 
varying magnetic field and tensor light shift used to drive 
the atomic spins. Our experience suggests that it may 



be possible to achieve far more precise control of the 
atomic ground state by driving the ensemble with radio- 
frequency and microwave radiation in the presence of a 
static bias magnetic field. This would eliminate opti- 
cal pumping as a source of decoherence and loss, fur- 
ther reduce the modest loss of state purity, and make 
the evolution almost completely reversible. It could also 
allow robust state estimation without invoking extrane- 
ous parameters and the positivity requirement for Pml, 
and this could shorten our data processing time from cur- 
rently ~ 1 / 2 hour to perhaps a few seconds on a desktop 
computer. At that point, further improvements in the al- 
gorithm and computing hardware — for example the use 
of programmable logic devices 15] — could plausibly al- 
low quantum state estimation and state-based feedback 
control in real time. 
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